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ABSTRACT 


This  report  contains  two  papers  on  the  dynamic  response  of  aircraft  operating  from  damaged  and  repaired  runways.  In 
the  first  paper  the  response  of  a  simplified  representation  of  an  aircraft  to  two  discrete  disturbances  is  analysed  to  see  how 
the  second  disturbance  modifies  system  behaviour  caused  by  the  first  disturbance.  The  second  paper  provides  a 
mathematical  model  which  can  be  used  for  calculation  of  the  dynamic  response  of  aircraft  structures  operating  on  rough 
surfaces;  a  comparison  is  made  between  theoretical  predictions  for  a  YF16  aircraft  and  typical  measurements  from 
frequency  response  tests. 


RESUME 


Deux  rapports  sont  inclus  dans  cette  publication  qui  traite  de  la  Reponse  Dynamique  d’un  avion  operant  sur  Piste 
Endommagee  et  Reparee.  Dans  le  premier  papier  est  analysee  la  reponse  a  une  representation  simplifiee  d’un  avion  ayant 
subi  deux  perturbations  quantifies  afin  de  voir  comment  la  seconde  perturbation  peut  modifier  la  conduite  du  systeme 
generee  par  la  premiere. 

Le  second  papier  traite  d'un  modele  mathematique  qui  peut  etre  utilise  pour  le  calcul  de  la  reponse  dynamique  des 
structures  d'un  avion  operant  sur  des  surfaces  reparees  mais  rugueuses;  la  comparaison  sera  faite  entre  les  predictions 
theoriques  pour  un  avion  YF 1 6  et  les  mesures  significatives  prises  a  partir  d'essais  de  reponses  en  mode  de  frequence. 
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INTERPRETATION  IN  TERMS  OF  THE  RESPONSE  OF  A  ONE  DEGREE-OF-FREEDOM  OSCILLATOR 

TO  TWO  SUCCESSIVE  DISTURBANCES 
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Air  Force  Wright  Aeronautical  Laboratories 
Wright- Patterson  Air  Force  Base,  Ohio  45433,  USA 

ABSTRACT 

This  paper  explains  the  dynamic  response  of  an  aircraft  that  taxies  over  two  arbitrary  disturbances,  under  the  assumption 
that  the  aircraft  can  be  represented  as  a  linear,  one  degree  of- freedom  system  .  That  analysis  produces  the  concept  of  the 
BUMP  MULTIPLIER  which  explicitly  and  simply  determines  whether  a  second  discrete  disturbance  will  amplify  or  attenuate 
the  response  from  a  first  disturbance.  The  BUMP  MULTIPLIER  also  simplifies  the  understanding  and  presentation  of  the 
results.  While  the  assumptions  are  very  severe,  the  resulting  formulas  can  be  very  useful  to  guide  more  elaborate  nonlinear 
calculations  or  to  plan  test  programs. 

1.  INTRODUCTION 

The  problem  of  aircraft  dynamic  response  to  taxiing  over  rough  surfaces  has  been  a  topic  of  research  for  many  years.  For 
the  most  part,  the  research  has  been  limited  to  predicting  and  measuring  the  dynamic  response  over  the  (nearly)  random 
roughness  imposed  by  the  terrain  or  by  wear  and  tear  on  runways  and  taxiways.  Within  the  last  several  years,  however, 
concerns  have  arisen  within  the  Defense  agencies  of  the  NATO  countries  about  the  safety  of  aircraft  operations  over  the 
discrete  disturbances  which  can  arise  from  bomb-damage  to  the  runways.  Those  concerns  also  extend  to  the  dynamic 
response  due  to  taxi  over  repaired  runways  and  repetitive  aircraft  operations  on  the  (potentially)  yielding  surfaces. 

As  a  result  of  those  concerns  the  United  States  Air  force  instituted  program  HAVE  BOUNCE  which  is  performing  taxi 
tests  over  simulated  (relatively  mild)  runway  damage  and  repairs  for  several  USAK  combat  and  transport  aircraft,  tne  project 
also  develops  computer  programs  to  predict  the  dynamic  response  to  the  simulated  runway  profiles  .  Other  NATO  nations 
are  performing  similar  test  and  analysis  programs  on  their  aircraft. 

HAVE  BOUNCE  considers  the  computer  programs  to  be  validated  when  they  produce  satisfactory  comparisons  with  the 
experimental  results  from  flight  (taxi)  tests.  Then  the  project  uses  the  validated  computer  programs  to  extrapolate  from 
the  relatively  mild  test  conditions  to  more  severe  operational  cases.  Since  the  computer  programs  account  for  the  nonlinear 
properties  of  the  landing  gears  and  tires,  the  extrapolation  of  their  results  beyond  their  validated  range  of  parameters  is 
always  open  to  some  question. 

Because  the  taxi  test,  programs  have  proven  to  be  very  expensive,  difficult  to  control  and  repeat  and  (sometimes)  danger¬ 
ous,  the  USAF  also  created  the  Aircraft  Ground  Induced  Loads  Excitation  (AGILE)  facility  which  measures  the  dynamic 
response  of  operational  aircraft  to  damaged  and  repaired  runways  within  the  controlled  conditions  of  the  laboratory.  The 
project  supports  an  operational  aircraft  on  its  tires  on  massive  hydraulic  shakers  and  drives  the  shakers  vertically  to  represent 
the  vertical  events  of  the  aircraft  taxiing  over  damaged  and  repaired  runways.  Earh  (of  the  three  integrated)  shakers  can 
sustain  a  static  weight  of  50,000  lb,  can  displace  amplitudes  of  10  in,  can  impose  dynamic  forces  up  to  50,000  lb  and  can 
be  driven  sinusoidally  (frequencies  up  to  25  Hz),  randomly  or  to  follow  prescribed  discrete  motions.  In  its  first  major  test, 
agreement  between  the  AGILE  tests  and  HAVE  BOUNCE  taxi  tests  for  an  operational  A-7D  aircraft  was  excellent. 

All  three  evaluation  methods  computer  programs,  HAVE  BOUNCE  taxi  tests  on  operational  aircraft  and  AGILE  tests 
on  operational  aircraft  have  been  dominated  by  one  major  consideration-the  nonlinearities  in  the  landing  gear.  As  a  result, 
nearly  all  of  the  computations  have  been  done  with  numerical  time- integration  of  the  nonlinear  differential  equations  of 
motion.  The  taxi  tests  and  AGILE  tests  also  have  been  forced  to  adopt  a  tedious  approach  of  repetitive,  trial  and  error  test 
cases,  again  because  the  aspect  of  nonlinearity  has  prevented  the  superposition  of  simple  disturbances  to  synthesize  more 
complex  responses. 

In  this  paper  we  contend  that  the  nonlinearities  do  indeed  strongly  influence  the  computational  and  test  results,  especially 
the  exact  levels  of  the  loads  obtained.  However,  the  qualitative  response  and  the  selection  of  speeds,  bump  heights,  and 
bump  sparings  which  produce  large  dynamic  responses  ought  to  be  predictable,  for  the  most  part,  by  simpler  linear  methods 
Nonlinear  calculations,  taxi  tests  and  AGILE  tests  all  ought  to  be  preceded  by  a  substantial  amount  of  linearized  calculations 
which  can  be  done  rapidly  and  can  yield  much  physical  insight  into  those  conditions  which  produce  extensive  dynamic 
response.  A  clever  analyst  may  Bbe  able  to  find  the  simplicity  and  intuitive  understanding  in  seemingly  complex  time 
histories,  which  in  fact  may  be  not  much  more  than  superpositions  of  many  relatively  simple  events. 

The  purpose  of  this  paper  is  to  review  those  linear  methods,  to  show  how  they  yield  an  understanding  of  complex  time 
histories  and  how  they  can  be  used  to  plan  nonlinear  calculations,  taxi  tests  and  AGILE  tests.  The  paper  illustrates  the 
principles  by  treating  the  response  of  a  linear  one  degree  of  freedom  oscilllator  as  it  taxies  over  two  successive  discrete 
disturbances,  introducing  the  concept  of  the  BUMP  MULTIPLIER. 
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2.  THE  RESPONSE  OF  A  ONE  DEGREE  OF  FREEDOM  OSCILLATOR  TO  A  SINGLE  DISTURBANCE 


Assume  a  one  degree  of  freedom  oscillator,  with  damping  less  than  the  critical  value,  receives  some  excitation  over  a 
period  of  time,  but  that  the  excitation  stops  at  time  t  =  tt. 


For  times  after  ()t  when  the  response  is  decaying  freely,  the  acceleration  response  can  be  written: 

5(t)!  {A|  sin  [ut(/ —  t| )]  +  fli  cos  [u»(t  —  <i)]}  (2.1) 

t>ti 


where: 

(i  =  the  time  the  excitation  ends 
w  =  damped  frequency 
iw  =  damping 

A\,B\  =  Constants  which  depend  on  a, ut,  the  excitation  and  the  initial  conditions 

Note  that  the  damping  parameter  o  above  is  not  quite  the  same  as  q,  the  frequently  used  fraction  of  critical  damping, 
which  comes  from  the  analysis  of  a  classical  one  degree-of  freedom  oscillator.  The  product  aui  controls  the  exponential  decay 
of  the  damped  system,  perhaps  as  observed  experimentally.  The  use  of  the  parameter  o  allows  us  to  refer  the  damping  to 
the  observed  damped  frequency  ut  rather  than  the  fictitious  undamped  frequency  u iq.  In  the  special  case  of  the  classical  one 
degree-of  freedom  oscillator  a  =  q / ^1~— "f* . 

The  decaying  acceleration  response  also  can  be  written  as: 

5(0|t>l|  =  R|e"""(<',|)sin|w(<  -<i)  +  d>,]  (2.2) 

where: 

Rl  =  yjA\  +  B * 
tand>,  =  B\/A\ 

We  loosely  refer  to  R \  as  the  potential  amplitude  of  the  acceleration  response.  It  is  an  upper  bound  on  the  amplitude  of  the 
acceleration  response  to  a  single  disturbance.  The  phase  shift  <t>\  depends  only  on  A \  and  Dt  and  will  therefore  be  different 
for  various  forms  of  the  excitation. 

For  small  damping  (o  <<  1)  the  behavior  of  the  acceleration  response  will  be  dominated  by  the  term  sin|u)(f  —  fi)  +  <£|j 
in  Equation  (2.2).  so  we  would  expert  its  local  maxima  and  minima  to  be  obtained  from  solutions  of 

w(t  -  ti)  +  0,  (2n  -  1)*;  n=  1,2,3,...  (2.3) 

However,  th  ■  term  e  au‘l  f|  I  causes  a  shift  in  the  values  of  time  for  which  the  local  maxima  and  minima  of  the  acceleration 
response  occur.  In  fact,  the  third  derivative  of  the  displacement  (derivative  of  the  acceleration)  is: 

rf^')-||>li  =*.«  ’’W1  +  o2ros[w(f  -  t,)  +  (0,  +  M)  (2.4) 

where: 

tan  t>;  =  a 

Therefore,  the  local  maxima  and  minima  of  the  decaying  acceleration  response  will  occur  at  the  values  of  time  for  which 

w(f  -  f,)  =  (2n  -  I)*  -(<*>,+  Ai ):  n  =  1,2,3,...  (2.5) 

The  additional  phase  shift  ,  will  be  small  for  values  of  damping  that  are  small  with  respect  to  the  critical  value,  (o  <<  1) 
Note,  however,  that  the  first  phase  shift  qi ,  need  not  be  small. 

We  can  use  these  results  to  obtain  an  even  better  upper  bound  on  the  amplitude  of  the  acceleration  response.  If  we  plug 
the  above  value  of  *■(/  -  f,)  into  Equation  (2.1)  for  c  ( # )  (  |  ,we  find  that  the  extreme  values  for  c(f)  (  (  equal: 

*!_«  "I* 20  11 a  (|>I  *  *'  •]  ■  n  =  1  2  3 
vAl  +  ft1 


I t>^tro«»  trnltm  nrnifs  for  ti  —  J . 
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S.  TWO  SUCCESSIVE  DISTURBANCES 


Now  suppose  the  single  degree  of  freedom  oscillator  receives  a  subsequent  excitation  over  another  period  of  time  and 
that  excitation  stops  at  time  I  =  l2. 


•  • 
2 


If  there  had  been  no  previous  response  from  the  first  disturbance  the  acceleration  response  to  the  second  disturbance 
would  have  been: 

5(<)'  ,  =  e  t3>  {A2sin[uj(t  -  <2)J  +  S2  cos  M<  -  <2)]} 

t  >t 2 

However,  because  of  the  presence  of  the  decaying  response  to  the  first  disturbance,  the  acceleration  response  to  the  combined 
disturbances  must  be  written: 


5(<)  ,  ,  =e  1,1  (A,  sin  [w(t  -  /,  )J  +  B,  cos  [w(<  -  <1))}  + 

r.  {.42sin[u/(f  -  t2)|  +  B2cos[u>(<  -  f2)]} 

For  convenience  in  manipulating  the  terms  in  Equation  (3,1)  we  abbreviate: 

e,  =  e  Qu,('  '•> 

S,  =  sin  [u»(<  -  t,)j 
C,  =  cos  |u;(f  -  (,)] 

Then  the  acceleration  response  to  the  combined  disturbances  is: 


5(0|(>ll  =  +B1C,)  +  e,(A,S,  +  B,CI)  (3.2) 

The  trick  is  write  Equations  (3.1)  and  (3.2)  with  respect  to  the  time  of  the  most  recent  disturbance,  t2.  To  that  end  we 
write: 

t-t,  =(<-tj)  +  (t2-ft) 

We  need  the  additional  abbreviations: 

e2i  =  e <i» 

52,  =  sin  |ut(<2  -  f , )] 

C2,  -  cos[w(f2  -  t,)] 

To  obtain: 

«t  =  e2e2) 

5,  =  S2C21  +  C2S2\ 

C|  =  C2C2(  -  S2S2i 

The  acceleration  response  to  the  combined  excitations  becomes 

2(/)|f>ej  =',ia;s2  +  b;c3)  (3.3) 

where: 

=  A2  +  e2t (^|C2i  -  B\S2X) 

B\  —  B2  +  e2|(MiS2i  +  B\C2x) 

Following  the  same  procedure  we  used  for  the  single  disturbance,  the  acceleration  response  to  the  combined  disturbances 
can  be  rewritten  as: 

5(0|l>tl  =  °w<t  (,)  M*  -  **)  +  <h]  (3-4) 

where. 

Bj  =  \Ja?  +  B27  (3.5) 

=  v'M|  +  B\)  +  2e2|  \C2X(AiA2  +  B,B2j  +  S2I(A,B2  -  BxAt)\  +  t\x(A]  +  B})) 


tanrii  -  +  eJl('4|‘‘’J)  +  BiC2i) 

a;  a2  +  c,,(A|C2I  -  R|-V||) 


1-4 


Equation  (S.S)  [or  R2,  the  potential  amplitude  oj  the  acceleration  response  to  the  combined  disturbances,  is  one  of  the 
major  findings  of  this  paper.  Much  of  the  subsequent  work  here  will  be  concerned  with  finding  the  conditions  which  maximize 
and  minimize  R2 

As  we  did  for  the  single  disturbance,  we  can  differentiate  Equation  (3.4)  with  respect  to  time  to  search  for  the  times  for 
the  local  maxima  and  minima  of  the  acceleration  response.  We  obtain: 

d^r\t>ti  -  +  a2cos[u;(t  -  t2)  +  (02  +  <52)]  (3.6) 

where 

tan  62  =  a 

W  e  see  the  same  phase  shift  ^  ^i)  in  the  times  for  local  local  n*axima  and  minima  of  the  decaying  acceleration  response 
to  the  combined  disturbances,  which  will  occur  when 

-•(' -'2)  =  (2n- 1)'  -(02  +  62):  "=1.2.3,...  (3.7) 

As  we  did  for  the  single  disturbance,  we  can  use  these  results  to  obtain  an  even  better  upper  bound  on  the  amplitude 
of  the  acceleration  response.  If  we  plug  the  above  value  of  jz(t  -  f2)  into  Equation  (3.4)  for  c(t)  (  (  ,  we  find  that  the  local 

extreme  values  for  c(f)  (  (  equal: 

»  °l<2n  Oj  ;n=l,2.3,... 


We  have  seen  how  to  find  the  times  for  local  maxima  and  minima  of  the  decaying  acceleration  response,  assuming  we 
know  .4j ,  ,  A2,  B2.  o,  .*•.  t !  and  t2.  However,  we  are  searching  for  the  best  and  worst  possible  runway  profiles,  so  the  most 

critical  aspect  is  to  find  the  values  of  (2,  the  time  of  the  second  disturbance,  which  will  locally  maximize  and  minimize  the 
potential  amplitude  R2.  We  differentiate  Equation  (3.5)  for  R2  with  respect  to  f2  and  set  the  result  to  zero  to  obtain: 

■Si  [ ( -4 1 A 2  +  B,B2)  +  q(A,B2  -  B|A2)]  + 

1 1  [o(A|.42  +  BxB2)  -  (A\B2  —  fl|A2)j  +  r»(A2  +  Bx)e2 1  =  0 


~Rn  sin  [w(f,  -  t , )  +  0„|  +  =  0  (3.9) 

vl  +  a2 

where 

.  _  a(A[A2  +  B\B2)  -  (AxB2  -  B\A2) 

(AxA2  -f  B\B2)  4-  o[A\B2  —  B]  A2) 

The  exact  solution  for  the  time  delays  (l2  -  f t )  which  locally  maximize  and  minimize  R2  would  require  a  numerical  or 
graphical  solution  of  Equations  (3.8)  or  (3.9),  however  for  small  damping  we  would  expert 

u(t2  -  tx)  nn  -  tl’l2\  n=  1,2,3...  (3.10) 

We  will  give  the  determination  of  the  time  delays  (t2  -  (,)  an  exact  treatment  in  the  next  Section. 

4  THE  DUMP  MULTIPLIER 

Recall  that  /?,  represented  the  potential  amplitude  of  the  decaying  acceleration  response  to  the  first  disturbance  and 
that  R2  represented  the  potential  amplitude  of  the  acceleration  response  to  the  combined  disturbances,  where  in  each  rase 
we  measured  time  from  the  time  of  the  most  recent  disturbance.  H  e  rail  the  ratio  R2/Rx  the  DVMP  MULTIPLIER,  since  it 
defines  the  extent  to  which  the  second  disturbance  amplifies  (or  attenuates)  the  response  to  the  first  disturbance  The  BI'MP 

MULTIPLIER  is 

R2  _  (.42  +  Of)  +  2e21  |f21(.4,.42  +  B  02)  +  N21(.4,ZJ2  -  B,.42|]  +  e2(,42  +  B‘)  ,  , 

.4f  +  ZJj*  ‘  (41) 

To  assist  in  the  interpretation  of  the  BUMP  MULTIPLIER  we  add  another  set  of  abbreviations: 


5  — 


1-5 


The  potential  amplitudes  and  phase  angles  for  the  responses  to  the  first  disturbance  and  the  combined  disturbances  become: 

=  |A,|Ul  +  f?  (4.2) 

tan  <p\  =  f  | 

1*2  —  I -42 1  ^  (1  +  (|)  +  2(  ^e2\  )  [^2l(l  +  <lf2)  +  -*>2 1  ( * 2  _  '  I  )]  +  (  ^  P21  )  (*  +  >U  I4-3) 

(2  +  (4*C2l)(^2l  +  1 1  C*2 1  ) 

tan  0>  = 

1  +  (  ^2  e2l  )  (C  21  -  <l*Sl) 

The  BUMP  MULTIPLIER  becomes: 


f?2  _  A2  (1  +  '2)  +  2(  A2e2l)  K'2I  ( 1  +  (1<2)  +  -S’2I  (*2  -  'l  )1  +  (  A^en)2(i  +  (f 

1  A,  |\ 


1  +  f  J 


(4-4) 


Note  in  Equation  (4.4)  that  the  magnitude  effects  are  contained  mostly  in  the  term  (}'  ,  whereas  the  spacing  effects  are 
in  the  radical.  The  spacing  effects  are  dominated  by  the  terms 

S2l  =  sin  [u>(<2  -  *i)] 

C2\  =  COs[u2(/2  -  <,)] 

and  those  terms  are  always  modulated  by  the  combination 

e.tl  -  e  <M<s  <l) 

A 2  21  -  A2 

When  we  use  these  abbreviations  in  Equation  (3.9)  to  find  the  time  delays  (l2  -  (|)  which  locally  maximize  and  minimize 
the  potential  amplitude  R2,  we  obtain 


!  ^2 

i*7l\ 


-^-|sin[w(<2  -fj)  +  4>l2\  +  7  "  y-  =  0 

t  +  c,  vl  +  o' 


(4.5) 


where 


tan0l2  =  ^1^V-Vr4-  Z<~\ 

(i  +  « 1  < 2 )  +  o(f 2  “  <l) 


A  first  approximation  for  o  <<  1  for  the  time  delays  would  be: 

—  #i)  =»  rz7T  —  0i2  ;n  =  I,2,3,... 

To  get  a  second  approximation  we  search  for  small  angles  dn  such  that 

^(<2  -  I,)  =  ntr  -  V>,j  -  f)n 

With  this  change  of  variables  Equation  (4.5)  takes  a  new  form 

np°0n 


C„  sin  = 


Vi  +  a2 


where 


C„  =  (-l)”^1"”  f| 

A\ 


\  1  + 1 


If  we  expand  Equation  (4.8)  for  small  d„,  keeping  only  first  order  terms,  we  obtain 

o 

I’n 


CnV 1  +  o2  -  a2 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4-10) 


This  then  is  the  (hoped  for)  small  correction  we  use  in  Equation  (4.7)  to  find  the  second  approximation  for  the  time  delays 
('2~'l). 


1 
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5.  USE  OF  THE  AVERAGE  SPEED 

We  have  made  no  assumption  of  a  constant  taxi  speed  between  the  two  disturbances;  in  fact  we  will  show  in  this  Section 
that  the  spacing  effects  are  dominated  by  the  average  speed  between  them. 

If  /  is  the  distance  between  the  two  disturbances,  the  average  speed  is 


V’  = 


l 

h-h 


Then  we  can  express  the  term 


-•('2  ->,)  =  7  =  A 


(5.1) 


(5.2) 


The  usual  terminology  for  A  —  ljj/V ,  based  on  the  instantaneous  speed,  is  the  reduced  frequency.  Therefore  A  -  Ij/V  is  the 
reduced  frequency  based  on  the  average  speed  between  the  two  disturbances. 


Amplitude  and  Du  trip  \t  ultiplier 

Equations  (4.2)  and  (4.3)  for  the  potential  amplitudes  and  phase  angles,  and  Equation  (44)  for  the  Bl'MP  Mt'LTIPLIER 
remain  unchanged,  except  that  now  we  note 

n‘\ 


e2i  =  « 


St 


sin  A 
cos  A 


Best /W  orst  Runway  Profile, 

The  process  for  finding  the  exact  and  approximate  solutions  for  the  time  delays  which  maximize  and  minimize  the  potential 
amplitude  Rj  of  the  decaying  acceleration  response  to  the  combined  disturbances  remains  unchanged,  except  Equation  (4.5) 
now  becomes: 

•■*2  1  +  ']  I  \  ,  , .  \  .  n  „  r>A  _  „ 

4  \  ,  >  si»l  t  A  +  t  I  2  )  4-  . f  —  l) 

-a,  \  i  +  '  ,2/  jm? 


(5.3) 


The  first  approximation  for  the  time  delays  becomes  an  equal  ion  for  A: 

A  *  nir  -  t  |2  ;  n  =  I.  2. 3.. 


(5.4) 


The  process  of  finding  exact  solutions  and  second  approximations  for  the  time  delays  also  remains  unchanged,  except 
that  Equation  (4.7)  becomes 

X  =  nit  -  vi2  -  (5.5) 


6  SPECIAL  C  ASE  SIMILAR  DISTURBANCES 

We  now  define  similar  disturbances  as  discrete  disturbances  that  have  the  same  shape  but  differ  only  in  magnitude 
and/or  sign.  Examples  would  be  the  entire  family  of  infinite  ramps  or  a  family  of  sine  waves  of  the  same  wavelength  but 
varying  heights.  The  a,,umption  of  similar  disturbances  i,  not  a  very  limiting  one.  In  fart,  nearly  all  of  the  profile s  tested 
in  the  HAVE  BOI’SCE  program  and  all  of  the  MATO/AdARD  profile,  ran  he  broken  down  into  sequenres  of  similar  ramp 
disturbances.  For  linear  systems  with  zero  initial  conditions  similar  disturbances  will  produce  similar  responses,  and  when 
disturbances  are  similar  f|  =  (2  -  < 

Amplitude 

Crider  the  assumption  of  similar  disturbances.  Equation  (4.2)  and  (4.3)  for  the  potential  amplitudes  and  phase  angles  become: 


R 1  7  i  A 1 ;  \  1  t  <2 
tan  O  |  =  r 
,A, 


Rj  .4/  \  1  *  1  4  2('i<V1||Ci,  ♦  !  ^r2i) 


(6.1) 


(6.2) 


tan  <i>i 


'  +  (  Ji,r21  )(-s2i  +  21 ) 


I  «■  (  Ij'-nHCz,  -  r.Sil 


D  a  mp  ,\f  till  i pi  ier 

Equation  (4.4)  for  the  Bl’MP  Mt’LTIPLIER  reduces  to; 


Ri 


■».  \ 


2(  ,  1  r2|  )(  21  +  (  ,  1  C2|  ) 
.4  2  ■ 1 2 


(6.3) 


▼ 


tWff  IT 


l  ■  |  ■  «  ■  ■  ■  . 


Figures  (6.1)  and  (6.2)  illustrate  the  BUMP  MULTIPLIER  versus  A  with  the  damping  as  a  parameter  f<*r  * t ^ •  *p«-  ^ 
case  of  equal  (or  opposite)  disturbances  at  constant  speeds. 

The  major  conclusion  to  hr  drawn  is  thr  dependence  of  the  Bl  SIP  Sfl'LTIPLIF.R  on  the  nondnnt  n.-otriai 
speed  or  frequency)  parameter 

a  -  -  M  -  I~iv 

The  spacing  of  the  disturbances,  the  average  speed  of  the  vehicle  and  the  natural  frequency  all  comhin*  info  ll-t  u**t  }h,r,,tt 
the  reduced  frequency  based  ott  the  average  speed.  Also,  while  the  Bt'StP  Sll'LTIPLIER  clearly  ilinumsli* »  with  ‘  f<>*  .«rj ;»  >'  f  .• 
values  of  n .  the  values  of  A  that  protluee  the  local  minima  and  maxima  art  seen  to  hr  fairly  n  *ak  funetotn.-  >>f  <> 


Best /Worst  Runway  Profiles 

In  order  to  find  the  time  delays  that  locally  maximize  and  minimize  the  potential  amplitude  R1  hquaMnn  ■>  3  r»du»** 


where 


The  definition  of  (  n  reduces  to 


•A, 


,(A 


02 


or 

v  1  ♦  fP 


tan  t'i 2  _  <> 


<n  ( 


<*  i  2  I 


At 

A{ 


i, 


4 
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Tables  (6.1)  ami  [6  2)  give  values  of  t he  first  approximation,  second  approximation  and  exact  solution  ft»r  *  wMrf;  io<a.l> 
maximize  arid  minimize  the  Bl’MP  MULTIPLIER  for  the  same  special  rase  of  equal  i>»r  opposite:  disturbances  and  ron*?an* 
speeds. 


The  major  conclusions  to  hr  drawn  are 

(a)  Again  the  sparing  that  maximizes  nr  minimizes  the  arcrlenition  rrsjMinsr  drprrids  irry  nearly  on  the  damping 
eter,  ct, 

(b) thc  First  approximation 

\  -  nr  -  i' a  ^  nr  n 

is  an  excellent  approximation  to  thr  exart  solution  for  reasonably  small  values  of  dwaping 


7  EXAMPLE  THE  CLASSICAL.  UNDER  DAMPED  SPRING  MASS  DAMPER  TAXIING  OVER  TWO 
RAMPS 

We  consider  the  example  of  a  classical,  single  degree  of  freedom  oscillator  that  encounters  two  ramp  disturbances  The 
disturbances  are  separated  by  a  distance  /  and  occur  at  times  /|  and  tj.  respectively  The  taxi  speeds  I'j.l’j  are  not  necessarily 
equal  at  the  time  of  the  encounters,  nor  are  the  ramp  angles  0j.0j. 


</(/)  -  vet  u[t) 

u(t)  unit  step  function 


1-8 


By  solving  the  ordinary  differential  equation  of  motion  for  the  displacement  in  response  to  the  ramp  inputs  (with  zero 
initial  conditions)  and  then  differentiating  those  results  twice  with  respect  to  time  we  find  the  various  parameters  to  use  in 
Equations  (2.1)  and  (3.1): 


A,  =  V,^0$,(  1  -2c2) 

B,  =  2Y,^O0,  cv  1  -  C2 

C  =  ratio  of  damping  to  critical  value,  c/2mui0 
u.'o  =  undamped  natural  frequency,  \j  k/m 
ui  =  damped  frequency,  woy  1  —  C2 

o  =  c/'/l  -  C2 

We  also  note  the  nondimensional  reduced  frequencies  : 

i  _ 

lu  / 

A  =  =  A0y  1  -  c2 


(7.2) 


(7.3) 


and  the  relationship: 


am  =  cu>o 


(7.4) 


Note  the  distinction  between  the  hypothetical  undamped  frequency  u>0  and  the  actual  damped  frequency  ui.  Note  also  the 
distinction  between  f,  the  fraction  of  critical  damping  as  related  to  the  hypothetical  undamped  frequency  uiq,  and  a  which 
we  use  to  relate  damping  to  the  actual  damped  frequency  u. 


so  that 


For  purposes  of  illustration  we  pick  the  fictitious  undamped  natural  frequency  to  be 

=  27t  rad/see 
fa  =  w0/2tr  =  1.0  Hz. 

and  we  pick  the  damping  value 

C  =  c/2mwo  =  0.1 

V'l  -  C2  =0.9950 
n  =  c/v'l  -  C2  =  0.1005 

Note  that  this  means  the  actual  damped  frequency  is: 

v  =  1.99tr  rad/ sec;  f  =  0.9950  Hz 


All  members  of  the  family  of  ramp  inputs  are  similar.  Therefore,  For  every  ramp  input  (regardless  of  speed  V’,  frequency 
jj  or  slope  0)  the  similarity  parameter  >  is: 


D  2cv/'l  -c2 

=  c,  =  <  =  .  =  -  _  ,  = 


2<r2 


tan  0|  =  0.2030 


(7.5) 


The  second  phase  shifts  for  the  location  of  the  local  maxima  and  minima  of  the  decaying  acceleration  response  will  be 
identical: 

6,  =  f2  =  S  =  sin  1  (  =  0.1002  rad  =  >.739n 
Figure  (7.1)  illustrates  the  two  phase  shifts  <t>\  and  (i  versus  the  damping  ratio 


i  Amplitude 

j 

i  The  potential  amplitudes  of  the  decaying  acceleration  responses  arc  obtained  by  applying  the  definitions  in  equations  (7.2) 

j  (7.4)  to  the  equations  for  similar  disturbances.  Equations  (6.1)  (6.5). 

i  Figure  (7.2)  illustrates  the  decaying  acceleration  response  of  the  classical  spring  mass  damper  for  a  typical  single  ramp 

input  at  a  constant  speed  of  V'  =  10.0  l/eer  ,  where  the  length  units  are  in  any  convenient,  consistent  system  Figure  (7.3) 

illustrates  the  sensitivity  to  speed  by  plotting  the  response  for  a  range  of  speeds  1'  =  10,11 . 20  //see.  Note  that  the 

amplitude  of  the  response  to  the  single  disturbance  grows  monotonically  with  increasing  speed.  Figure  (7.4)  illustrates  a 
curious  feature  in  the  dependence  of  the  response  on  damping  for 

C  =0.1, 0.2,. ...0.5 

at  a  constant  speed  of  V  -  10.0  l/.ier.  Large  values  of  damping  actually  increase  the  peak  acceleration  response  at  the 
earliest,  instants  of  time.  These  effects  arr  the  result  of  the  term 

1  D,  =  2V',u.'ofl,  Cv^l  -  f2 


t 


J 


m 

$ 


& 

'‘Vi 


tjj 

ra3 


which  reaches  its  maximum  value  at  c  -  'j2 


Figure  (7.5)  illustrates  the  decaying  acceleration  response  for  traversing  two  equal  (but  opposite)  ramps,  separated  by  a 
distance  of  20.0  /.at  a  constant  speed  of  10.0 //see  Note  that  the  maximum  amplitude  of  the  acceleration  response  to  the 
combined  disturbances  at  a  speed  of  10.0 //sec  is  actually  smaller  than  the  response  to  the  first  disturbance  alone  at  the  same 
speed.  Figures  (7.6)  and  (7.7)  show  that  this  will  not  always  be  the  case  By  changing  the  constant  speeds  to  l'  =  8.02  and 
13.46  l/acc  for  the  same  ramp  geometry,  we  see  that  the  response  to  the  combined  disturbances  can  be  markedly  greater 
than  the  response  to  the  first  disturbance  alone  Figure  (7.8)  further  illustrates  the  sensitivity  to  speed  by  plotting  the 

response  for  a  range  of  speeds  l’  =  10,  11 . 20  / /sec.  While  the  amplitude  of  the  response  to  the  single  disturbance  grew 

monotonically  with  increasing  speed,  the  amplitude  of  the  response  to  the  combined  disturbances  displays  a  much  more 
complicated  structure. 

Figure  (7  9)  shows  the  damping  effect  on  the  acceleration  response  to  the  combined  disturbances  Note  how  the  maxi- 
mums.  zeroes  and  minimums  occur  at  nearly  the  same  periodic  values  of  A,  regardless  of  the  value  of  the  critical  damping 
ratio  f . 

Bump \f  ul/ipj  irr 

Reverting  to  the  general  case  of  non-constant  speeds  and  ramp  angles,  the  potential  amplitude  for  the  dynamic  acceleration 
response  to  the  combined  disturbances  is  given  by 

e>  _  t-  .  I*  .  .  ^  •>  I  \  .  c-  ^  /  1  itfirH  17  Si 


v'r-ol*2i\  l  +  2  (  )  cj.r’z,  +  (  [.£ 


Cj,  =  t  “-t'j  'll  -  e  “•»  =  f  0  IOOSA 
(  -  co -  /,)]  rr  cos  A 

The  BI'MP  M1LTIPLIER  from  Equation  (6  3)  is 


Ri  '  i9i  .  ,  ,  /  v  i#ieJi  \  /  1  i^i^zi  \ 

ft,  “  \  1  ‘  1  VJ9l  ) ( 11  l  r,9,  / 


The  term  in  Equation  17  7)  gives  the  pure  magnitude  effect  of  the  two  disturbances  The  radical  gives  the  sparing 

effect  The  dominant  term  is  Cj,  =  cos  A.  which  is  modified  by 

I  1  1*1  \  _  _  /  *  1*1  1  _  „i  _ .  /  *  1*1  \  _  0  I00S4 


Equation*  (7  6)  ami  (7.7)  art  very  powerful  results  uhtrh  relate  the.  potential  amplitude  Rj  and  the  Bl'\1P  MI  L- 
TIPLIER  to  the  truit  an  tan  eon .«  speed*  the  arerage  speed  l*.  the  ramp  angle  *  the  damping  parameter  n  and  the 

average  reduced  frequency  A  l+i'V 

Best /to  or.it  Runu  ay  Profile a 

Recall  that  the  equation*  which  determine  the  best /worst  runway  profiles  or  speeds  wrere: 

"1*  'in  I  V  .  i-tJ)  .  0  (6.4) 

4  l  v  1  *  " 

where 

tan  V|  j  -  o 

For  tin*  example 

*in  1  v  0  1002  rad  739“ 

For  constant  speeds  and  equal  (hut  opposite)  ramp*  Equation  16  4)  reduces  to 

*'  '  ♦  *in  (  A  ♦  t,  w )  0  (7  8) 
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a  First  Approximation 

The  first  approximation  to  tin-  time  delay*  that  locally  maxmii/c  anil  minimize  the  acceleration  response  to  the  combined 
disturbances  is: 

i  >  nr  -  t'lj:  n  =  1,2,3...  ( o .4 ) 

or 

A  -  3  04 14.6  1830. 9.3246.  12.4662, . . . 

b.  Exact  Solution 

In  obtaining  the  exact  solution  for  the  tune  delays,  the  values  of  for  n  -  1.2.3...  are: 

<  „  -  2 1  (44.  480  9:  11270;  310.7(H) 

The  corresponding  exact  (rapidly  diminishing)  values  of  ,fn  are 

.(„  4  779  •  10  '.  2  064  .  10  '.8  920  •  10  h.  3.855  -  10  7  .  .  . 


j  and  the  exact  values  of  \  are 

i  A  3  0366.6  1832.9  3246.  12  4662, .  . 

These  small  values  of  („  illustrate  why  the  first  approximation  is  such  an  excellent  approximation. 

■  c  Second  Approximation 

5  The  approximate  solutions  fur  small  •  would  be 

L  4  n 

vlt  o2  -  n‘ 


(4.10) 


These  approximations  to  f„  are 


I 


4  774  ■  10  -2.064  •  10  4. 


They  differ  from  the  exact  values  only  for  n  1,  and  then  very  slightly  In  addition,  both  the  exact  solutions  and  the  second 
approximations  are  very  small  corrections  to  the  first  approximation: 


A  = 


«ir  -  ti’u  =:  ntr  -  a 


(5.4) 


Earlier  in  this  section,  when  we  found  the  speeds  8.02 and  13.46 //.sec  that  increased  the  response  to  two  equal  (but 
opposite)  disturbances  separated  by  a  distance  20.0/.  we  used  the  results  in  Table  (7.2).  Table  (7.2)  shows  that,  for  a  %  0.1, 
the  exact  solutions  for  the  values  of  the  reduced  frequency  A  =  l+’/V  that  maximize  and  minimize  the  acceleration  response 
to  the  combined  disturbances  -re: 


Maximize:  A/it  =  0.9447.2.9558.4.9616.. 
Minimize:  A/ir  =  1.9853.3  9774.5  9731. 


For  a  length  /  =  20.0  /.  and  noting  that  /  =  j./2ir  =  0  995//;,  these  results  translate  to  the  following  speeds: 

Maximize:  V  =  42. 13,  13.46.  8.02  //see. ..  . 

Minimize:  V  =  20.05.  10.01.6,66 //.sec. .  .  . 

Note  that  Table  (7.1)  or  (7.2)  requires  the  use  of  the  damped  reduced  frequency.  A  =  /ui/V'  rather  than  the  the  fictitious 
undamped  reduced  frequency  A0  =  l^n/\' . 


8  APPLICATION  TO  NONLINEAR  CALCULATIONS  AND  TEST  PROGRAMS 

a.  Threr  Pri rin plr* 

The  first  set  of  basic  ideas  to  keep  in  mind  when  when  using  these  results  to  plan  nonlinear  calculations,  HAVE  BOUNCE 
(taxi)  tests  or  AGILE  tests  is  that  a  useful  building  block  is  the  infinite  ramp,  that  two  infinite  (opposite)  ramps  can  combine 
to  produce  an  AUARD  Dump,  and  that  two  (opposite)  AGARD  bumps  can  combine  to  produce  an  AGARD  Repair  Mat. 


Ill 


Second,  recall  that  for  two  disturbances  separated  by  a  distance  /,  the  best/worst  combinations  tend  to  occur  when: 

-  Iw  2j r// 

A  =  =  — - —  ss  nx  —  a  ;n  =  1,2,3... 

V  V 


V  ~  2'4)  ;n  =  1,2,3... 


(8.1) 


Third,  the  potential  amplitudes  of  the  acceleration  response  of  a  classical  one  degree  of  freedom  oscillator  to  a  single 
infinite  ramp  and  to  two  combined  ramps  are: 

/?,  =  vlWte|#,| 


Rt  =  ^o\f[v2e2y  +  2iylel)(v2e2)e2lc2l  +  (v,«,)2 


where: 


-V^1+2(^1)cm  +  (^) 

e21  =  e-°"(tj-ti)  =  e-«A 
C2i  =  cos  [w(t2  -  ti)j  =  cos  A 


(7.6) 


b.The  Response  on  the  Initial  Slope  oj  a  Repair  Mat 

Beginning  with  the  basic  building  block  of  the  infinite  ramp  we  have  shown  that  the  acceleration  response  is  given  by: 

S(0it>li  =  ft,e  -““(‘-‘•)«in|w(t-l,)+*,|  (2.2) 


The  initial  acceleration  response  at  the  first  corner  will  be  given  by: 

i\t=t+  =  R|sind>i 

We  have  also  shown  that  the  local  extreme  values  of  the  acceleration  response  are  given  by: 

-■A_ ;n  =  1,2,3,... 

vTTS’ 

and  that  the  peak  values  occur  when 

w(fpeo*  _  <i)  =  (2n  -  1)  j  -  (<t>\  +  6|)  ;n  =  1,2,3,... 
This  translates  to  distances  where  the  peaks  occur  of: 


V  f  „• 

ip ««*  -  *i  =  -  |(2n  -  l)j  -  (^i  +  6|)l  ;n  =  1,2,3,... 


_ (2n-l)V 

Xpeak  x\  ~ 


2(<t>\  +  6 

A  7n  . 


ir(2n  -  1) 


;  n  =  1,2,3,... 


(8.2) 


(8.3) 


where  V  is  the  average  speed  over  the  initial  slope. 

These  are  the  expressions  for  the  peak  amplitudes  of  the  acceleration  response  and  the  locations  of  those  peak  responses 
on  the  initial  slope  of  an  AGARD  Repair  Mat.  Conversely,  if  the  length  of  the  initial  ramp  is  ia ,  we  know  that  the  peaks  of 
the  acceleration  response  will  occur  somewhere  on  the  ramp  for  speeds: 

V  < - r~  ;n  =  1.2.3,...  (8.4) 

<*"  -*)[•- 

Otherwise  the  peak  acceleration  responses  on  the  initial  slope  -.v ill  occur  at  the  second  corner.  At  that  point  the  value 
becomes: 

5,  =R,e 


sin  (  'y  +  d>,  ) 


=  Rxe 


(8.5) 
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c.  The  Flat  Area  o f  a  Repair  M at 

We  assume  that  angles  0,,02  are  equal  and  opposite  and  that  l\ ,  the  length  of  the  ramp,  is  small  enough  so  that  the 
speed  over  the  initial  slope  is  constant,  \\  =  Vj.  Then  the  acceleration  response  to  the  combined  disturbances  is: 


where 


-  R2'  a‘,(<  sin  [u>(t  -  t2)  +  ^2] 


(3.4) 


f?2  -  | v' 1  -  2e2iC2i  +  e^, 


tan  <t>  2  = 


(  -  e2i (521  +  (C2{) 

i  -  e2i(C2i  -  <521) 


The  acceleration  response  just  after  the  second  corner  is: 


h  ,  R2sin<t>2 

t  f  9 


The  local  peak  values  of  the  acceleration  response  on  the  flat  part  of  the  repair  mat  are  given  by: 

123 

y/TT^ 


(8.6) 


and  the  peak  values  occur  when 


l\,» 


This  translates  to  distances  where  the  peaks  occur  of: 

[ (2n  -  l)x 


xptak  x2 


(2 n  -  1)V 

xprak  x2  ~ 


—  (<t>  2  +  £2) 
2(<t>2  +  <§2) 


jr(2n  -  1) 


;  n  =  1, 2, 3, . . . 


;  n  =  1,2,3,... 


(8.7) 


where  in  this  case  V  is  the  average  speed  over  the  flat  part  of  the  repair  mat.  If  the  length  of  the  flat  part  is  /2,  we  know 
that  peak  acceleration  responses  will  occur  somewhere  on  the  flat  part  for  speeds: 


T/  ^  _  4/^2  _  _  .  n  n 


(8.8) 


d.  Obtainingthe  Infinite  RampData  f  rom  the  Test  Results  for  an  AGARD  BUMP 
Because  of  the  impossibility  of  experimentally  developing  an  Infinite  Ramp ,  it  will  be  more  practical  to  excite  the  oscillator 
with  an  AGARD  Dump  and  then  infer  what  the  response  would  be  to  an  Infinite  Ramp.  We  assume  that  we  have  excited 
the  oscillator  with  an  AGARD  Bump  and,  therefore,  we  will  know  e2i,  52i  andC2 \  and  will  have  measured  R2,a,ui,<t>2  in  the 
equation: 

S(<)jf>(j  =  R2e  sin  [u>(t  -  t2)  +  <£2j  (3.4) 

The  potential  amplitude  R |  ran  be  obtained  from: 


Then  the  phase  lag  4>\  can  be  obtained  from: 


R 1  = 


5, 

2C21C21  +  e*] 


C2,5 
1  -  e2lC 


tan  <t>2  +  Q 
1  -  Q  tan  <t>2 


tan  <t>  1  =  f 


(8.9) 


(8.10) 

(8.11) 
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e.  Guidelines  for  Nonlinear  Calculations,  AGILE  Tests  andTaxi  Tests  We  begin  by  calculating  or  measuring  the  ac 
celeration  response  to  an  AGARD  Bump 


over  a  range  of  speeds  V  and  angles  9.  Since  the  response  will  undoubtedly  not  be  purely  in  a  single  degree  of  freedom,  we 
must  process  the  test  data  to  obtain  separate  values  of  It 2.  o,  and  <S>2  for  each  degree  of  freedom.  Then  for  each  speed,  angle 
and  degree  of  freedom  we  calculate: 

e21  =  e  C2]  =  cosj^ft,  -  / , )j 

and  we  use  Equations  (8.9)  (8.11)  to  calculate  R\  and  ip ,  A  good  test  of  our  assumed  linearity  is  to  form  The  values 
should  be  approximately  the  same  for  each  degree  of  freedom,  regardless  of  speed  V  or  angle  $ 

The  next  step  should  be  to  test  the  linear  result  that  the  Best/Worst  AGARD  Bump*  will  be  those  for  which 

A  -  =  ntr-o;n  =  l.2.3. 

Since  1 1  is  fixed  by  the  AOARD  geometry,  we  ran  accomplish  this  variation  by  choosing  the  speeds  to  be 


We  ran  interpret  the  final  slope  of  the  AGARD  Repair  Mai 


as  just  the  negative  of  the  initial  slope;  with  the  only  distinction  that  it  occurs  at  a  distance  of  I,  +  Ij  after  the  initial  slope 
Then  we  ran  search  for  the  Best/Worst  length  of  the  repair  mat  by  setting 

;  ('1  +  /»)-• 

A  -  -  n*  -  o;  n  1 , 2.  3. 

V 

In  this  rase  we  have  both  the  average  speed  V  and  the  length  l2  to  use  as  variables 

Now  we  note  that  the  total  length  of  the  AGARD  Repair  \fat  is  21]  +■  l2  and  assume  that  another  repair  mat  is  placed  1 
distance  1,  behind  the  first  mat 

Therefore,  to  look  for  the  Best/Worst  sparings  we  set 

,  .2/,  -  1 2  4  /,).• 


I 


9.  CONCLUSIONS 


We  have  treated  the  dynamic  response  of  an  aircraft  taxiing  over  runway  disturbances,  under  the  assumption  that  the 
gross  aspects  of  the  dynamic  response  can  be  found  in  the  analysis  of  a  linear,  one  degree  of  freedom  system,  excited  by  two 
successive  disturbances.  We  have  found: 

a.  There  is  a  great  deal  that  can  be  learned  about  the  governing  physics  for  aircraft  dynamic  response  to  taxi  over 
damaged  and  repaired  runways  by  examining  the  results  of  calculations  with  relatively  simple,  linear  models. 

b.  The  seemingly  complicated  time  histories  can  be  merely  superpositions  of  relatively  simple,  time-phased  events. 


r  Relatively  simple  expressions  are  available  for  the  potential  amplitude  (an  upper  bound)  of  the  acceleration  response 
excited  by  one  or  two  disturbances.  In  the  (not  too)  special  case  of  similar,  disturbances  separated  by  a  distance  /,  with 
nonconstant  speeds  ami  ramp  angles,  the  expression  for  the  potential  amplitude  R 2  is: 


where 


Rl  -  l  i^il  ^  1  + 


,  I  U,0te2,  \ 

1 1  V202  I 


Cn  + 


I  l'i ^1  e2i  \2 

\  r202  ) 


(”21  =■  cos  A 
A  u- fi¬ 


ll.  The  effects  of  disturbance  spacing  and  variable  taxi  speed  are  controlled  by  the  reduced  frequency,  based  on  the 
average  speed  between  disturbances: 


e  One  need  not  actually  calculate  the  time  histories  to  find  the  best/worst  profiles  and  speeds,  but  ran  use  the  expressions 
for  the  potential  amplitude  Rt  and  the  DUMP  MULTIPLIER  ^ 

f  To  maximi/e/mimmize  dynamic  response  a  good  approximation  for  A  is 

A  =  titr  -  ft 


g  While  damping  obviously  controls  the  dynamic  response  to  the  disturbances,  the  critical  speeds  and  disturbance 
-parings  are  weak  functions  of  damping 

These  results  ran  ra-ilv  be  extended  from  two  disturbances  to  an  arbitrary  number  of  disturbances  and  multiple 
degree  of  freedom  system-  with  multiple  landing  gear 

1  The  results  of  calculations  based  on  these  linear  methods  should  be  compared  with  results  from  flight(taxi)  tests, 
AGILE  test-  and  nonlinear  calculations  Tin-  i-  not  to  say  that  the  linear  results  should  be  relied  on  to  predict  detailed 
loads:  rather  the  question  should  be  tin  ihr  simple  /incur  models  prrdirl  Ihr  rriliral  speeds  and  spanngs  so  that  we  ran  use 
them  to  <jut tie  our  ft  jtt  program#  anti  nonlinear  #011111011# 
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Table  7.1  The  Disturbance  Spacing  \/ir  to  Maximize  and  Minimize  the  Response  to 


Two  Equal  Disturbances 
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Table  7.2  The  Disturbance  Spacing  A/x  to  Maximize  and  Minimize  the  Response  to 
Two  Equal  and  Opposite  Disturbances 
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AN  EXPERIMENTAL- ANALYTICAL  ROUTINE  FOR  THE  DYNAMIC  QUALIFICATION  OF 
AIRCRAFT  OPERATING  ON  ROUGH  RUNWAY  SURFACES 
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SUMMARY 

A  mathematical  model  to  be  used  as  a  basis  for  analytical  investigations  to  predict  the  dynamic  struc¬ 
tural  response  of  flexible  aircraft  operating  on  rough  runway  surfaces  is  presented.  It  is  shown  how  the 
structural  parameters  included  in  the  aircraft  generalized  equations  of  motion  are  determined  in  a  ground 
vibration  test  on  the  real  aircraft  structure  and  in  additional  tests  on  components  of  the  undercarriage. 
The  validation  of  the  developed  mathematical  model  is  achieved  by  a  comparison  of  typical  results  from  fre 
quency  response  tests  and  calculations  performed  on  a  YF-16  prototype  fighter  aircraft.  Finally,  the  way 
in  which  the  developed  mathematical  model  can  be  used  in  combination  with  various  systematic  test  pro¬ 
cedures  for  the  dynamic  qualification  process  of  aircraft  operation  on  damaged  /  repaired  runways  is  in¬ 
dicated  . 
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circular  eigenfrequency  of  the  r-th  eigenmode 


value  related  to  r-th  global  aircraft  eigenmode  (r  =  1,2,. ...n) 

value  related  to  the  ji-th  landing  gear  degree  of  freedom  (u  =  1,2 . m) 

value  related  to  a  tire  parameter 

value  related  to  a  landing  gear  parameter 

value  related  to  a  suspension  system  parameter 

value  related  to  a  landing  gear  strut 

denote  first  and  second  derivatives  with  regard  to  time 


1.  INTRODUCTION 

For  many  years,  the  NATO  countries  have  been  dealing  with  the  problem  of  evaluating  the  capatibilities 
of  their  aircraft  to  operate  from  rough  runway  surfaces.  This  problem  is  a  primary  anticipation  in  modern 
warfare  tactics,  which  foresee  an  early  destruction  of  the  opponent's  airfields. 

To  guarantee  the  readiness  of  aircraft  even  in  the  case  of  damaged  runways,  engineers  have  con¬ 
centrated  on  the  task  of  correctly  predicting  the  dynamic  structural  response  of  aircraft  when  taxiing  over 
rough  runway  surfaces  and  to  establish,  based  on  these  data,  guidelines  (criteria)  for  the  (rapid)  repair 
of  damaged  runway  surfaces. 

In  order  to  determine  realistic  limits 
of  aircraft  operational  capabilities  on  the 
ground,  extensive  experimental  air¬ 
craft  taxi  investigations  <  Figure  1 )  have 
been  performed  [1.2).  These  tests,  al¬ 
beit  associated  with  many  difficulties, 
have  proven  able  to  realistically  deter 
mine  the  aircraft  taxi  capabilities  on 
rough  surfaces  But  since  this  testing 
is  also  very  costly,  taxi  testing  will  at 
ways  be  restricted  to  a  minimum. 

On  the  other  hand,  the  many  ana 
lytical  investigations  performed  in  the 
past  indicate  that  .  for  calculation  of  re¬ 
alistic  dynamic  response  data,  it  is  nec 
essary  to  consider  in  the  mathematical 
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model  of  the  aircraft  structure,  both  the  aircraft  rigid  body  and  flexible  eigenmodes  as  well  as  the  non¬ 
linear  and  frequency-dependent  characteristics  of  the  undercarriage  [3].  This  requirement  is  easily  for¬ 
mulated,  but  is  more  difficult  to  realize.  For  the  engineer  involved  in  structural  dynamics,  numerous 
questions  result  therefrom,  such  as  how  the  aircraft  structural  model  should  be  built  in  order  to  allow 
easy  handling  of  the  structural  nonlinearities,  and  how  the  (modal)  parameters  of  the  flexible  aircraft 
structure  as  well  as  the  characteristics  of  the  undercarriage  can  be  realistically  determined. 

Closer  insight  into  the  overall  problem  indicates  the  need  to  develop  a  practical  mixed  experimental- 
analytical  routine  which  allows  experimental  determination  of  realistic  structural  parameters  consistent  with 
a  well-adapted  mathematical  structural  model,  to  be  used  as  a  basis  in  the  dynamic  response  calculations. 

This  report  gives  a  description  of  a  newly  developed  experimental  -  analytical  method  which  can  be 
adopted  as  a  routine  in  the  qualification  process  of  aircraft  for  its  operation  on  damaged /repaired  runways. 
The  practical  applicability  of  the  method  is  demonstrated  by  the  investigations  performed  on  a  prototype 
YF-16  fighter  aircraft.  In  the  following  chapters  and  paragraphs,  first  a  model  of  the  aircraft  structure  is 
presented.  Moreover  it  is  shown  how  the  (structural)  parameters  included  in  the  aircraft  equations  of  mo¬ 
tion  can  be  experimentally  determined  in  various  tests  performed  on  the  aircraft  structure  and  on  its  un¬ 
dercarriage.  Finally  the  developed  mathematical  structural  model  will  be  validated  by  comparing  calculated 
and  measured  frequency  response  results  of  typical  structural  data. 

The  work  described  in  this  paper  was  performed  by  the  author  at  the  Air  Force  Flight  Dynamics  Lab¬ 
oratory  of  the  Wright-Patterson  Air  Force  Base  in  Dayton  (Ohio)  during  a  one-year  stay  as  "visiting  sci¬ 
entist".  The  author  wishes  to  acknowledge  the  authorities  of  the  Structures  and  Dynamics  Division  for 
having  provided  the  opportunities  to  investigate  the  developed  routine  on  a  fighter  aircraft  of  the  US  Air 
Force.  Gratitude  is  also  due  to  the  many  people  at  the  Flight  Dynamics  Laboratory  who  graciously  pro¬ 
vided  their  support  during  the  various  tests. 


3.  THE  AIRCRAFT  STRUCTURAL  MODEL 

The  generalized  equations  of  motion  of 
flexible  aircraft,  as  depicted  in  Figure  2  , 
were  derived  from  the  Lagrange  Equations 
[4],  which  are  expressed  in  the  form 
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for  a  set  of  n  generalized  coordinates  qr> 
related  to  n  generalized  degrees  of  free¬ 
dom. 

Considering  n  generalized  degrees  of 
freedom  of  the  flexible  aircraft  structure 
and  m  additional  generalized  degrees  of 
freedom  for  the  translatory  motion  of  the 
undercarriage,  the  generalized  equations  of 
motion  of  the  entire  aircraft  can  be  formu¬ 
lated  as  follows: 


Figure  2:  Physical  structural  aircraft  model 
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I 

with 

(3)  AxL(t)  =  X0  '  Aq(t) 

as  the  landing  gear  deformations , 

I  <4>  Vl)  =  Z0L’P(t) 

1  i 

as  the  runway  roughness  encountered  by  the  landing  gears ,  and  j 

(5)  *(t)  =  ZQ  •  p(t) 

as  the  (rigid-body)  displacements  at  discrete  points  of  the  structure  induced  as  the  structure  encounters 
runway  roughness.  At  any  given  time,  the  total  (absolute)  displacement  at  a  well-defined  point  P  of  the 
aircraft  structure  is  given  by  superposition  of  the  elastic  and  rigid-body  mode  response  and  of  the  forced 
1  displacement  due  to  the  runway  roughness.  Thus  we  can  write 

(6)  xp(t)  =  ♦p-q(t)  +  Zop  •  p(t)  ,  (P  =  1,2 . N)  . 

As  indicated  by  Equation  (2),  consideration  of  the  landing  gear  degrees  of  freedom  as  "additional"  de¬ 
grees  of  freedom  entails  that  the  structural  parameters  related  to  the  overall  elastic  aircraft  structure  are 
completely  separated  from  the  landing  gear  parameters.  This  formulation  of  the  equations  of  motion  is  ad-  ' 

vantageous  if  the  dynamic  structural  parameters  of  the  landing  gears  have  to  be  modified  often  during  cal¬ 
culation,  as  for  example  when  nonlinear  landing  gear  characteristics  must  be  considered. 

In  the  following  chapters  the  way  in  which  the  various  structural  parameters  included  in  the  gener¬ 
alized  equations  of  motion  can  be  experimentally  determined  will  be  shown. 


4.  GROUND  VIBRATION  TEST 

The  aim  of  a  ground  vibration  test  is  to  provide  the  modal  structural  data  of  a  flexible  structure. 
These  parameters  are 

•  the  eigenfrequencies  f  , 

•  the  normal  mode  shapes  $r  (eigenmodes)  , 

•  the  generalized  masses  Mrf  , 

■  the  generalized  (modal)  damping  coefficients  Drr- 

Common  ground  vibration  test  procedures  such  as  the  classical  "Phase  Resonance  Method",  require  the 
structure  to  be  nearly  linear.  This  is  definitely  not  true  for  an  aircraft  structure  supported  by  its  un¬ 
dercarriage.  Investigating  this  aircraft  structure  by  a  ground  vibration  test  would  entail  that  the  meas¬ 
ured  modal  data  are  a  function  of  the  amplitude  level  of  the  external  excitation  applied  to  the  aircraft 
structure  in  the  test.  As  a  consequence,  experimental  determination  of  a  set  of  normal  mode  shapes  is 
not  possible  for  this  aircraft  configuration. 

But ,  if  the  basic  aircraft  configuration  is  modified  such  that  the  relative  piston  motion  in  the  landing 
gear  struts  is  fuUy  suppressed,  i.e.  that  the  landing  gear  struts  are  locked,  the  aircraft  structure  can  be 
regarded  as  being  nearly  linear.  The 
equations  of  motion  of  the  aircraft  con¬ 
figuration  with  locked  landing  gear 
struts  can  be  derived  from  Equation  (2) 
by  setting  Aq  =  0  .  Moreover  it  is  as¬ 
sumed  that ,  for  simulation  of  the  cor¬ 
rect  boundary  conditions  in  the  ground 
vibration  test ,  the  aircraft  is  supported 
at  its  nose  and  main  landing  gear  tires 
by  a  soft  (pneumatic)  spring  system 

( cc  «  cT  ,  u  =  1,2 . m)  with  very 

low  damping  properties  (dgy  -*  0).  The 
YF-16  ground  vibration  test  setup  [5] 
is  depicted  in  Figures  3  and  4  . 


The  condition  Cg^  <<  can  be  easily 

satisfied  when  inflating  the  aircraft  tires 
to  a  "high"  pressure.  With  Q(t)  as  the 
vector  of  the  generalized  external  excita¬ 
tion  forces,  the  equations  of  motion  of 
the  aircraft  (ground  vibration  test)  con 
figuration  with  locked  landing  gears  can 
be  formulated  as  follows: 

(7)  M  <j(t)  ♦  [[K  ♦  +1  eg  #L]  * 

♦  i  [O  ♦  «IT  ♦L]]q(t)  =  Q(t)  . 

Application  of  the  classical  phase  reso¬ 
nance  ground  vibration  test  method  [6] 
to  this  system  allows  the  eigenfrequencies 
f  ,  the  generalized  masses  M  (Table  1) 
and  the  normal  mode  shapes  (Figures 
5  and  6)  to  be  determined  directly. 


MODE 

EIGEN  FREQUENCY 
fr  [Hz] 

GENERALIZED  MASS 
Mrr  [K?  m2] 

DAMPING  COEFFICIENT 

*r  t%DcritJ 

RB  1 

Rigid  Body  Pitch 

0.621 

1890.0 

/ 

RB  2 

Rigid  Body  Heave 

1.410 

4611.0 

/ 

RB  3 

Rigid  Body  Roll 

2.019 

106.0 

/ 

S  1 

Symm.  Bending 

4.62 

870.0 

1.6 

S  2 

Symm.  Missile  Pitch 

7.01 

275.0 

1.5 

S3 

Fuselage  Vert.  Bending 

12. 19 

317.0 

1.7 

A  1 

Asymm.  Missile  Pitch 

5.73 

1186.0 

1.7 

A  2 

Asymm.  Wing  Bending 

8.21 

346.0 

1.4 

A3 

Fin  Bending 

14. 10 

33.2 

1.9 

A  4 

Fuselage  Lat.  Bending 

1677 

110.0 

1.7 

Table  1:  Modal  parameters  of  YF-16  aircraft 


The  values  of  the  modal  damping  coefficients  are  fixed  or  derived  as  follows: 


for  all  rigid-body  modes 


and 
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for  the  flexible  modes,  0r  being  the  experimentally  determined  modal  damping  coefficients,  as  denoted  in 
Table  1.  This  approximation  of  the  structural  damping  is  permissible,  since  the  values  of  the  coefficients 
in  are  small  relative  to  flexible  eigenmodes . 

The  coefficients  of  the  (nondiagonal)  stiffness  matrix  K  are  determined  from 


(8) 


['“>rMrrJ  ‘  *1  CS 


which  is  possible  when  the  stiffnesses  Cg^  of  the  aircraft  suspension  system  are  known.  Determination  of 
these  stiffnesses  can  easily  be  achieved  in  a  static  loading  test  (Table  2). 

Thus,  after  completion  of  the  ground  vibration  test,  the  following  parameters  are  available: 

•  the  diagonal  generalized  mass  matrix  M, 

•  the  diagonal  modal  damping  matrix  D, 


Figure  6:  First  symmetric  wing  bending  mode 


Nose  LG*1 

Left  Main  LG*5 

Right  Main  LG  5 

cgl  =  1.34  •  105  N/m 

cg2  =  1.56  105  N/m 

cg3  =  1.50  105  N/m 

Landing  gear 

Table  2:  Measured  stiffnesses  of  the  YF-16  ground  vibration  test  suspension  system 


•  the  nondiagonal  generalized  stiffness  matrix  K, 

•  the  modal  matrix  4  and  its  submatrix  > 

•  the  diagonal  matrix  eg  containing  the  suspension  stiffnesses. 

A  review  of  the  equations  of  motion  of  aircraft  with  unlocked  landing  gears  (Equation  (2))  indicates  that 
all  parameters  related  to  the  overall  flexible  aircraft  structure  can  be  determined  from  the  ground  vibra¬ 
tion  test  investigations  as  described.  The  remaining  unknown  (nonlinear)  parameters  are  all  related  to  the 
undercarriage.  Experimental  determination  of  these  parameters  in  additional  tests,  performed  on  the  land¬ 
ing  gear  struts  and  tires  will  be  demonstrated  in  the  following  two  chapters. 


5.  LANDING  GEAR  STRUT  TEST 

A  major  problem  encountered  when  performing  dynamic  aircraft  taxi  response  calculations  is  how  to 
obtain  realistic  data  for  approximation  of  the  spring/damper  characteristics  of  the  landing  gear  (strut)  el¬ 
ements.  Usually  the  required  data  are  not  available,  since  the  investigations  performed  by  the  landing 
gear  manufacturer  normally  concentrate  on  (transient)  landing  gear  drop  tests  for  simulation  of  the  land¬ 
ing  impact. 
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Figure  7:  YF-16  main  landing  gear  strut 

test  setup 


For  the  performance  of  dynamic  response  calculations, 
however,  we  need  to  know  the  dynamic  stiffness  and 
damping  characteristics  of  the  landing  gear  strut  elements. 
These  characteristics  can  be  determined  by  (harmonically) 
cycling  the  landing  gear  strut  in  a  testing  machine  ( Fig¬ 
ure  7)  at  various  frequencies  and  amplitudes  of  its  dis¬ 
placement.  To  obtain  correct  landing  gear  data  for  various 
aircraft  gross  weights,  the  landing  gear  strut  tests  must 
be  performed  at  different  levels  of  the  static  strut  preload 
force . 

For  a  strut  preload  force  Fq  =  75  kN  ,  consistent  with 
the  static  loading  of  the  strut  installed  in  the  aircraft  in 
the  case  of  the  YF-16  configuration,  cycling  tests  were 
performed  on  a  two-stage  main  landing  gear  strut.  The 
experimentally  determined  values  of  the  (equivalent)  dy¬ 
namic  stiffness  [7]  are  depicted  in  Figure  8.  The  plotted 
curves  indicate  a  strong  nonlinear  but  a  nearly  frequency- 
independent  behavior  of  the  dynamic  strut  stiffness.  Fig¬ 
ure  9  depicts,  at  a  frequency  of  1  Hz  ,  the  equivalent  stiff¬ 
ness  Cg^,  of  the  strut  as  a  function  of  the  amplitude  of 
vibration  x  . 


Figure  8:  Experimentally  determined  nonlinear 
stiffness  behavior  of  the  YF-16  main 
landing  gear  strut 


Figure  9:  Experimentally  determined  equivalent 
dynamic  stiffness  of  the  YF-16  main 
landing  gear  strut 


The  dynamic  damping  characteristics  of  the  strut  were  obtained  from  measured  hysteresis  curves  [8], 
some  of  which  are  depicted  in  Figure  10.  A  graph  of  the  viscous  damping  factor  dg^,  is  shown  in  Figure 
11  ,  indicating  a  strong  nonlinear  and  frequency-dependent  damping  behavior  of  the  strut. 


The  thus  determined  strut  characteristics  c„_  and  dc_  (u  =  1 . m)  are  not  identical  to  the 

landing  gear  characteristics  Cj  ^  and  dj  considered  in  Equation  (2).  But  their  respective  values  are  in¬ 
terrelated  by  a  factor  or  a  nonlinear  function  g^  =  g  (AX|(j),  wh'ch  can  be  derived  from  the  landing 


( u  =  1 . m )  . 


gear  geometry. 

The 

following  interrelations  exist 

(9a» 

ci 

Ui 

CSTu  ' 

( geometry )  , 

(9b) 

d. 

Lm 

dSTu 

( geometry)  , 

Thus,  after  performance  of  the  landing  gear  strut  tests,  the  following  data  arc  available: 


Figure  10:  Measured  hysteresis  curves  on  the  YF-16  main  landing  gear  strut  at  a  static  preload  of  75  kN 
a)  f  =  0.25  Hz  b)  f  =  1.0  Hz  c)  f  =  3.0  Hz  d)  f  =  5.0  Hz 


the  landing  gear  stiffness  characteristics 
cLu  =  cLm(F0m’  f- AxLu)  •  (w  =  1 . m) 


•  the  landing  gear  (viscous)  damping  characteristics 

d.  =  d.  ( Fn  ,  f,  Ax.  )  ,  ( u  =  1 . m)  . 

Lu  Lm  0m  Lu 

Discrete  values  of  these  different  characteristics  are  the  coef¬ 
ficients  of  the  diagonal  matrices  C L  and  ,  respectively. 
With  regard  to  evaluating  the  coefficients  of  the  landing  gear 


mass  matrix 


it  has  to  be  mentioned  that  their  determina¬ 


tion  is  of  no  difficulty,  since  the  values  of  the  unsprung 
masses  m^^  (u  =  1 . m)  are  exactly  known  by  the  manu¬ 

facturer. 


6.  TIRE  TEST 

To  obtain  realistic  tire  stiffness  and  damping  data,  tests 
must  be  performed  on  spinning  tires.  Dynamic  tire  stiffness 
and  damping  data  can  be  experimentally  determined  in  a  test 
setup,  as  depicted  in  Figure  12.  The  operation  of  the  test 
setup  is  as  follows: 

The  runway  surface  is  simulated  by  the  rolling  surface  of 
a  wheel,  spinning  at  a  high  angular  velocity  and  driving  the 
aircraft  tire.  A  hydraulic  shaker  (harmonically)  cycles  the  tire 


|  12  mm 
0  ! 


Preload  =  75  kH 


2  3  <  5  Hz  6 

f - - 


Figure  11:  Experimentally  determined 
viscous  damping  factor  of 
the  YF-16  main  landing 
gear  strut 
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on  the  wheel.  Tests  are  performed  at 
various  wheel  speeds  V  =  ■  p  ,  various 

static  tire  preloads  F^  and  for  various 
frequencies  f  and  amplitudes  x  of  the 
dynamic  shaker  excitation.  Measuring  the 
dynamic  force  F  and  the  corresponding 
displacement  x  ,  allows  determination  of 
the  tire  stiffness  characteristic 

cT  =  F/x  =  c^,  (  Fq  ,  V  ,  f ,  x ) 

measurement  of  hysteresis  curves  allows 
determination  of  the  tire  structural  damp¬ 
ing  characteristics 

Figure  12:  Setup  for  dynamic  tire  testing  dT  d^,  (  Fq  .  V  ,  f .  x) 

Discrete  values  of  the  cT^  and  d^,^  characteristics  (p  =  l,...,m),  determined  for  the  different  tires  of 
the  nose  and  main  landing  gear  are  entered  as  coefficients  into  the  diagonal  matrices  e^,  and  .  re¬ 
spectively.  Thus,  after  completion  of  the  tire  tests,  the  matrices  e^,  and  d^  are  well  defined. 

Within  the  scope  of  the  investigations  performed  on  the  YF-16  aircraft,  tests  were  only  carried  out  on 
the  non-spinning  tires  of  the  nose  and  main  landing  gears.  The  investigations  concentrated  on  the  deter¬ 
mination  of 

•  the  (nonlinear)  load/stroke  tire  characteristics  from  which  the  (nonlinear)  stiffness  character¬ 
istic  can  be  derived , 

•  (quasi-static)  tire  hysteresis  curves  for  evaluation  of  the  structural  tire  damping. 

For  performance  of  the  tests,  the  different  tires  were  squeezed  in  a  testing  machine  (Figure  13).  As 
a  typical  result,  Figure  14  depicts  the  static  load/stroke  curve  of  the  main  landing  gear  tire.  It  can  be 
noticed  that,  at  higher  preloads,  the  tire  stiffness  is  nearly  a  constant.  Tire  nonlinearities  are  of  second¬ 
ary  importance  when  compared  to  the  highly  nonlinear  behavior  of  the  landing  gear  struts.  But,  should 
the  tire  load /stroke  curve  indicate  larger  nonlinearities,  the  methods  described  in  [7]  can  be  used  for  de- 


Figure  13:  Machine  used  for  the  static  and  Figure  14:  Measured  load/stroke  curve  of  a 

quasi-static  tire  testing  YF-16  main  landing  gear  tire 
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termination  of  an  equivalent  amplitude-dependent  stiffness 
Figure  15  shows  hysteresis  curves 


Cij  ( F  q  i  x )  . 


measured  around  a  fixed  static  tire  preload  FQ  =  30  kN . 
From  the  hysteresis  curves,  an  average  structural  damp¬ 
ing  value  of  5%  D  ^  was  determined. 


7.  DYNAMIC  RESPONSE  INVESTIGATIONS 

As  all  structural  parameters  denoted  at  the  left-hand 
side  of  Equation  (2)  are  known  after  completion  of  the 
ground  vibration,  landing  gear  strut  and  tire  tests,  dyn¬ 
amic  response  calculations  can  be  performed  for  a  well- 
defined  runway  roughness  input,  provided  the  matrix  at 
the  right-hand  side  of  the  equation  system  is  known.  At 
first  it  seems  that  the  discrete  (finite  element)  mass  ma¬ 
trix  of  the  aircraft  structure  must  be  known  for  determi- 

*L  Z0L 


Figure  15:  Measured  hysteresis  curves  on 
the  YF-16  main  landing  gear 
tire  at  a  static  preload  of  30  kN 

].  This  "problem"  can  be  avoided  when  writing 


nation  of  the  right-hand  submatrix  [t  m 
the  matrix  of  displacements  Tq  in  a  series,  as  a  superposition  of  the  rigid-body  mode  displacements 
determined  in  the  ground  vibration  test,  as  follows: 


UO)  *0  =  <»r-a  • 

the  coefficients  of  matrix  A  being  the  factors  of  the  linear  combination.  Considering  the  existing  ortho¬ 
gonality  relations  between  the  (measured)  eigenmodes ,  we  can  write: 

(11)  l*rmZQ  ♦  OLmLZ0L]  =  [$]  *  A  , 

M  being  the  diagonal  generalized  mass  matrix,  as  defined  at  the  left-hand  side  of  Equation  (2).  The  co¬ 
efficients  of  matrix  M  have  been  determined  in  the  ground  vibration  test.  Equation  (11)  indicates  that 
there  is  no  direct  external  excitation  in  the  elastic  aircraft  eigenmodes  by  the  runway  roughness.  The 
elastic  modes  are  all  excited  by  their  respective  rigid-body/elastic  mode  coupling  terms  at  the  left-hand 
side  of  Equation  (2). 


7 . 1  Dynamic  Response  Testing  with  AGILE 


By  virtue  of  the  Aircraft  Ground  Induced  Loads  Excitation  (AGILE)  shaker  test  facility  (Figure  16) 
now  available  at  the  Air  Force  Flight  Dynamics  Laboratory  at  Wright  -  Patterson  Air  Force  Base  in  Dayton 
(Ohio),  extensive  dynamic  response  testing  on  real  aircraft  structures  has  become  possible  [9].  This  fa¬ 
cility  allows  performance  of  special  investigations  which  could  not  be  realized  in  the  past.  For  instance  it 

is  possible  to  harmonically  excite 
the  aircraft  via  its  tires  up  to  frequen¬ 
cies  in  *  -  vicinity  of  30  Hz  and  to  per¬ 
form  tt^  s  in  a  locked  landing  gear  air¬ 
craft  configuration.  These  possibilities 
first  seem  to  be  of  no  importance  since 
neither  this  type  of  excitation  nor  this 
configuration  are  representative.  With 
regard  to  the  qualification  of  the  aircraft 
structural  model,  however,  this  special 
type  of  testing  is  advantageous. 

First,  excitation  of  the  structure  by 
a  (swept)  sine  forcing  function  allows  ex¬ 
citation  of  all  eigenmodes  in  the  frequen¬ 
cy  range  of  interest.  A  comparison  be¬ 
tween  measured  and  calculated  frequency 
Figure  16:  The  AGILE  shaker  test  facility  response  data  (transfer  functions)  clearly 


indicates  which  of  the  modes  are  correctly  modelled  and  which  are  not.  secondly ,  the  ability  to  test  the 
aircraft  in  its  locked  landing  gear  configuration  allows  verification  of  whether  the  flexible  aircraft  struc¬ 
ture  itself  -  without  consideration  of  the  dynamically  very  complicated  landing  gear  -  has  been  modelled 
correctly.  This  intermediate  verification  step  is  important  in  a  response  analysis,  since  there  definitely  has 
to  be  agreement  between  measured  and  calculated  data  for  the  aircraft  locked  landing  gear  configuration 
before  concentrating  on  the  far  more  problematic  (real)  configuration  with  unlocked  landing  gear  struts. 

A  survey  of  the  experimental  and  analytical  investigations  performed  on  the  YF-16  aircraft  in  both  of 
the  configurations  with  locked  and  unlocked  landing  gear  struts  is  given  in  the  next  two  sections. 

7.2  Locked  Landing  Gear  Configuration 
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A  full  description  of  the  AGILE  fre¬ 
quency  response  tests  performed  on  the 
YF-16  aircraft  (Figure  17)  is  given  in 
[10]  .  Aircraft  response  data  were  rec¬ 
orded  and  analyzed  for  a  series  of  dif¬ 
ferent  shaker  excitation  configurations, 
e.g.  excitation  of  the  aircraft  at  the  nose 
or  main  landing  gear  tires  only  or  simul¬ 
taneous  excitation  at  both  of  the  landing 
gear  systems.  The  aircraft  was  excited 
by  the  shakers  with  a  constant  amplitude 
swept  sine  signal  in  the  frequency 
range  from  0.6  to  15  Hz.  Plots  of  typi¬ 
cal  dynamic  response  data,  resulting  from 

a  test  with  simultaneous  excitation  of  the  Figure  17:  AGILE  testing  on  the  YF-16  aircraft 

nose  and  main  landing  gear  tires,  are 

plotted  as  dashed  curves  in  Figures  18  to  21.  Moreover  these  plots  depict  the  corresponding  results 
from  dynamic  response  calculations  based  on  Equation  (2)  with  consideration  of  Aq  E  0.  Only  symmetric 
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Figure  18:  Main  landing  gear  tire  deformation 
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Acceleration  at  the  fuselage  nose 


Figure  21:  Acceleration  at  the  left  wing  tip 
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aircraft  eigenmodes  were  taken  into  account  in  the  calculations,  i.e.  the  modes  RBI,  RB2,  SI,  S2  and  S3 
of  Table  1.  All  calculations  were  performed  on  the  basis  of  the  true  measured  structural  data,  with  no  up¬ 
dating  of  any  of  the  structural  parameters  with  regard  to  an  improvement  of  the  calculated  results. 

The  figures  indicate  good  overall  agreement  between  measured  and  calculated  data.  Especially  the  fre¬ 
quencies  of  the  modes  are  well  determined  by  the  calculation.  Some  discrepancies  between  the  amplitude 
data  are  due  to  the  "bad"  resolution  (0.15  Hz  interval  between  frequency  steps)  of  the  digital  Fast  Fou¬ 
rier  Transform  Analyzer  used  in  the  data  reduction  process.  The  broader  "peaks"  measured  in  some  re¬ 
sponses  (e.g.  in  the  12.5  Hz  mode  of  Figure  20)  are  probably  due  to  a  nonlinear  behavior  of  the  aircraft 
structure  which  was  excited  at  real  high  g- levels  during  the  response  tests. 

In  general,  the  agreement  between  measured  and  calculated  response  data  was  found  to  be  accurate 
enough  to  start  the  dynamic  response  calculations  on  the  aircraft  configuration  with  unlocked  landing  gear 
struts. 


7.3  Configuration  with  Unlocked  Landing  Gear  Struts 


Analogous  to  the  procedure  described  in  Section  7.2,  experimental  and  analytical  frequency  response 
investigations  were  performed  on  the  YF-16  aircraft  in  its  (normal)  unlocked  landing  gear  configuration. 
Figures  22  to  25  depict  typical  frequency  response  data  for  the  same  shaker  excitation  configuration  as 
investigated  in  Section  7.2,  but  for  a  higher  amplitude  level  of  the  excitation.  In  the  calculations  the  ex¬ 
perimentally  determined  nonlinear  landing  gear  characteristics  were  considered. 
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Figure  24:  Main  landing  gear  strut  deformation 


Figure  25:  Acceleration  at  the  left  wing  tip 


The  figures  show  that  there  is  reasonable  agreement  between  calculated  and  measured  results,  except 
in  case  of  the  acceleration  value,  depicted  in  Figure  25,  the  correlation  cannot  be  regarded  as  satisfactory. 
The  discrepancy  between  measured  and  calculated  curves  in  this  plot  is  mainly  due  to  the  fact  that  two 
asymmetric  modes  (at  6  and  8.5  Hz)  are  excited  by  the  symmetric  shaker  excitation.  While  performing 
subsequent  tests  on  the  landing  gears,  it  was  noticed  that  this  effect  has  to  be  attributed  to  the  different 
transfer  behavior  of  the  left  and  right  main  landing  gear  struts. 


8.  CONCLUSION 


The  many  investigations  described  in  this  report  have  shown  that  the  dynamic  (frequency)  structural 
response  of  aircraft  to  an  external  excitation  at  its  landing  gear  tires  can  be  correctly  calculated.  To  ob¬ 
tain  satisfactory  results,  however,  it  is  required  that  the  mathematical  model  of  the  aircraft  structure  be 
carefully  adapted  to  a  series  of  tests  allowing  experimental  determination  of  all  structural  parameters  in¬ 
cluded  in  the  generalized  equations  of  motion.  It  was  noted  that  especially  the  nonlinear-,  frequency-  and 
preload -dependent  characteristics  of  the  undercarriage  have  to  be  carefully  considered  in  the  calculations. 


Despite  the  extensive  investigations  already  performed,  further  effort  is  required  for  the  developed 
mathematical  model  to  be  used  in  time  domain  analyses  allowing  the  prediction  of  the  aircraft  structural 
response  to  u  (discrete)  runway  roughness,  e.g.  bumps,  repairs.  But  the  entire  philosophy  of  aircraft 
testing,  as  described  in  this  paper,  in  combination  with  the  developed  mathematical  structural  model  can 
be  used  as  a  basis  for  a  routine  allowing  dynamic  qualification  of  aircraft  for  operation  on  damaged /repaired 
runways.  Figure  28  depicts  in  a  block  diagram  the  entire  qualification  process  consisting  of  a  series  of 
systematic  analytical  and  test  investigations.  The  routine  makes  ample  use  of  the  AGILE  shaker  system  as 
a  qualification  test  setup.  For  the  validation  of  the  mathematical  structural  model  it  is  considered  to  be  of 
importance  to  perform  AGILE  tests  on  both  of  the  aircraft  configurations  with  locked  and  unlocked  landing 
gear  struts.  But  since  AGILE  cannot  simulate  spinning  tire  effects,  this  facility  does  not  have  the  capa¬ 
bility  to  fully  replace  the  taxi  field  tests.  Taxi  field  tests  on  real  aircraft  will  remain  necessary  to  com¬ 
pletely  demonstrate  the  aircraft  operational  capabilities  on  the  ground. 


Figure  26:  Block  diagram  of  the  aircraft  qualification  routine 
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14.  Abstract 

This  report  contains  two  papers  on  the  dynamic  response  of  aircraft  operating  from  damaged  and 
repaired  runways.  In  the  first  paper  the  response  of  a  simplified  representation  of  an  aircraft  to  two 
discrete  disturbances  is  analysed  to  see  how  the  second  disturbance  modifies  system  behaviour 
caused  by  the  first  disturbance.  The  second  paper  provides  a  mathematical  model  which  can  be 
used  for  calculation  of  the  dynamic  response  of  aircraft  structures  operating  on  rough  surfaces,  a 
comparison  is  made  between  theoretical  predictions  for  a  YETb  aircraft  and  typical  measurements 
from  frequency  response  tests 


